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ABSTRACT 

We investigate the non-spherical density structure of dark halos of the dwarf spheroidal (dSph) 
galaxies in the Milky Way and Andromeda galaxies based on revised axisymmetric mass models from 
our previous work. The models we adopt here fully take into account velocity anisotropy of tracer 
stars confined within a flattened dark halo. Applying our models to the available kinematic data of 
the 12 bright dSphs, we find that these galaxies associate with, in general, elongated dark halos, even 
considering the effect of this velocity anisotropy of stars. We also find that the best-fit parameters, 
especially for the shapes of dark halos and velocity anisotropy, are susceptible to both the availability 
of velocity data in the outer regions and the effect of the lack of sample stars in each spatial bin. 

Thus, to obtain more realistic limits on dark halo structures, we require photometric and kinematic 
data over much larger areas in the dSphs than previously explored. The results obtained from the 
currently available data suggest that the shapes of dark halos in the dSphs are more elongated than 
those of ACDM subhalos. This mismatch needs to be solved by theory including baryon components 
and the associated feedback to dark halos as well as by further observational limits in larger areas of 
dSphs. It is also found that more diffuse dark halos may have undergone consecutive star formation 
history, thereby implying that dark-halo structure plays an important role in star formation activity. 

Subject headings: galaxies: dwarf spheroidal galaxies - galaxies: kinematics and dynamics - galaxies: 
structure - Local Group - dark matter 


1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies in the Milky Way 
(MW) and Andromeda (M31) galaxies are excellent lab¬ 
oratories to shed light on the nature of dark matter 
because these galaxies are the most dark matter dom- 
inate d systems, with total dynamical-light ratios of 10 - 
1000 (lMateoll998HGilmore et alJl2007t IWolf et al.ll2010D . 
DSph galaxies also have the advantage that individual 
member stars can be resolved, thus, it is possible to mea¬ 
sure very accurate line-of-sight velocities for their mem¬ 
ber stars. Therefore, using these high-quality data, we 
are able to constrain their internal structure of dark ha¬ 
los in light of the currently standard A-donrinated Cold 
Dark Matter (ACDM) models. In particular, the stud¬ 
ies of dSphs are important to understand several contro¬ 
versial issues that ACDM models hold on galactic and 
sub-galactic scales. 

One of the outstanding issues in ACDM models is 
the so-called “core-cusp” problem: ACDM-based N- 
body simulations have predicted stron gly cusped pro¬ 


files i n the center of dark halos jNarorro^^I 


1997t iFukushige fc Makinol Il997t iDiemand et al 


1996 


2008; 


Ishivama et al.112013IL whereas the dark halos in the 


observed galaxies, especially dSphs and low surface 
brightness (LSB) galaxies, suggest a cored or sh al¬ 
lower cusped density profiles rtMoor 
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Recent studies of this problem have argued the possibility 
that a fraction of the observed dSp hs and LSBs can ac¬ 
tually have cusped density profiles (lHavashi et al.l 120041: 
Strigari et al.ll20l7i iBattaglia et alJl2011b iStrigari et all 


201411 . We note that some specific model assumptions 


(such as spherical symmetry) and degeneracies (such as 
between velocity anisotropy and mass distribution) yield 
uncertainties in the model results, thus it is yet unclear 
if the most or all of the dSphs indeed have cored central 
density profiles in their dark halos. 

Another example of the tension between ACDM pre¬ 
dictions and observations is the “too-big-to-fail” prob¬ 
lem that masses of dark halos associated in the dwarf 
satellites of the MW are significantly lower than those 
of the most massive subha los i n an MW-sized halo in 
the ACDM universe dBovlan-Kolchin et al.|[20lTl . 120121 ). 
This problem may be regarded as a rewriting of the 
“core-cusp” problem, where the presence of a core in 
the center of a halo tends to lower i ts mean mass den¬ 
sity inside the luminous par ts (e.g., iBrook fe Di Cintiol 
12014 lOgiva fe Burkertll2015l) . Thus, to solve both prob¬ 
lems simultaneously, a possible solution may rely on a 
transformation mechanism from a cusped to cored cen¬ 
tral density, e.g., through the effects of baryonic process. 
In recent years, high-resolution cosmological TV-body and 
hydrodynamical simulations have attempted to alleviate 
the above issues; inner dark halo profiles at dwarf-galaxy 
scales can be made cored due to energy feedback from 
gas outflows driven by star-formation activity of galax¬ 
ies such as radiation energy from ma ssive stars, stel¬ 
lar w i nds, and superno v a exp l osions dGovernato et all 
2012; jPenarrubia^et^L l2012t TGarrison-Kimme l et al.l 


2013: iBrooks fe Zolotov |2014; iDi Cintio et al.l 2014 


Madau et al.l 12014 lOgiva fe Moril 12014 lOnorbe et al.1 
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TABLE 1 

The Observational Data set for MW and M31 dSph Satellites 


Object 

No. of Stars 

M, 

(10 6 M 0 ) 

M v 

Dq 

(kpc) 

Dialf 

(pc) 

q' 

(axial ratio) 

T0.7 a 

(Gyr) 

Refeferences b 

MW dSphs 

Carina 

776 

0.38 

-9.1 ±0.5 

106 ±6 

241 ± 23 

0.67 ±0.05 

4.5 

1,3,7 

Fornax 

2523 

20.0 

-13.4 ±0.3 

147 ± 12 

668 ± 34 

0.70 ±0.01 

4.6 

1,3,7 

Sculptor 

1360 

2.3 

— 11.1 ± 0.5 

86 ±6 

260 ± 39 

0.68 ± 0.03 

12.7 

1,3,7 

Sextans 

445 

0.44 

-9.3 ±0.5 

86 ±4 

682 ±117 

0.65 ±0.05 

No data 

1,3,7 

Draco 

185 

0.29 

-8.8 ±0.3 

76 ±6 

196 ± 12 

0.69 ±0.02 

11.5 

2,7 

Leo I 

328 

5.5 

-12.0 ±0.3 

254 ± 15 

246 ± 19 

0.79 ± 0.03 

2.7 

2,3,7 

Leo II 

200 

0.74 

-9.8 ±0.3 

233 ± 14 

151 ± 17 

0.87 ±0.05 

6.8 

1,7 

M31 dSphs 

And I 

51 

3.9 

-11.7 ±0.1 

745 ± 24 

670 ± 30 

0.78 ± 0.04 

7.6 

5,6,7 

And II 

488 

7.6 

-12.4 ±0.2 

652 ± 18 

1230 ± 20 

0.90 ±0.02 

6.2 

4,6,7 

And III 

62 

0.83 

-10.0 ±0.3 

748 ± 24 

400 ± 30 

0.48 ±0.02 

8.8 

5,6,7 

And V 

85 

0.39 

-9.1 ±0.2 

773 ± 28 

350 ± 20 

0.82 ±0.05 

10.0 

5,6,7 

And VII 

136 

9.5 

-12.6 ±0.3 

762 ± 35 

770 ± 20 

0.87 ±0.04 

12.8 

5,6,7 

a This value 
fromlWeisz et 

is the lookback 
all I 20 TIL 

time at achieving 70% of 

current stellar 

mass of dSphs, 

and is estimated by available data taken 


b Refe rences: ( I )ljxrdlL-A_LLu zidimil l ion (1995); (2) Martin ct al. ( 2 0081 ) : ('3'I IWalker et all d2009cT l: Ml lHo et al.H2012f ): (5') IToIIerud et aLl 
(2012); (6) McConnachic & Irwin (2006); (7) McConnachic (2012) 


120151) . For instance, using the results from cosmolog¬ 
ical ACDM simulation of isolated dwarf galaxies that 
can reproduce observed stellar mass-metallicity relations, 
metallicity distribution func tions, and s t ar-for mation 
histories in the Local Group, iMadau et al.l (12014D inves¬ 
tigated whether both the core-cusp and too-big-to-fail 
problems can be mitigated based on supernova feedback. 
They found that in the simulated dwarfs, the total en¬ 
ergy injected by SNe II exceeds the minimum energy 
required for the core-cusp transformation at all times. 
Thus, they concluded that simulated dwarf galaxies can 
form a cored dark halo and have a maximum circular ve¬ 
locity that is lower than the case for dissipationless dark 
matter only, because supernova feedback can effectively 
reduce dark matter density in the inner part of its dark 
halo. In the light of suggestions from these simulations, 
dark halo structures can be related with the evolution of 
dwarf galaxies such as star-formation history (SFH) and 
chemical enrichment. 

The shapes of dark halos, especially those of subha¬ 
los associated with a host dark halo, provide important 
information on the dynamical evolution of subhalos, as 
well as the structure of stellar systems formed within 
their subhalos. This is because the shapes of subha¬ 
los are generally subject to the str ength and f req uency 
of tidal effects from the host halo dKuhlcn ct al.l 120071) 
and facilitate appropriate investigations about the va¬ 
lidity of some of the assumptions in th e mass modelin g 
of ste llar kinematics in dwarf galaxies (|Vera-Ciro et al.l 
l20ll) . IV-body simulations based on ACDM theory 
have res olved that dark halos h a ve generally triaxial 
shapes ( iJing fe Siitol 1200(11 120021 : lAllgood et, al.l 120061 : 
lHavashi et al.l I2007T) . Also, all of recent 7V-body simu¬ 
lations have predicted that dark subhalos, which are at 
the low-mass end of a mass function f or dark ha los as well 
as galaxies, are not str ongly triaxial dKuhlen et al.ll2007t 
I Schneider et al.l 1201 21) but rather statistica lly oblate, 
axisymmetric shapes (IVera-Ciro et alJ [2014). Photo¬ 
metric observations of dSphs in the MW and M31 
also show that their lig ht distributions on the sky 
are actually no n-spherical (|IFwin fe Hatzidimitrioulll995l : 

iMcConnachie fe Irwinl 120061) . Therefore, to describe the 


internal structures of non-spherical dark halos and stellar 
systems in dSphs, we should use non-spherical mass mod- 
el s for these galaxie s rathe r than spherical mass models. 

lHavashi fo Chibal (120121 hereafter HC12) constructed 
axisymmetric mass models, where each of luminous and 
dark matter distribution has a non-unity axial ratio, and 
applied these models to line-of-sight velocity dispersion 
profiles of six bright dSphs in the MW to constrain the 
non-spherical structure of their dark halos. They con¬ 
cluded that the shapes of these dark halos are very flat¬ 
tened (axial ratio of dark halo: Q ~ 0.3 — 0.4) for most 
of the sample dSphs. We note that their axisymmetric 
mass models assume the distribution function of stars in 
the form f(E,L z ) : where E and L z denote the binding 
energy and angular-momentum component toward the 
symmetry axis, respectively. Under this simplified as¬ 
sumption, velocity dispersions of stars, v\ and v%, in 
the R and z directions in cylindrical coordinates, re¬ 
spectively, are identical, or in other words, a velocity 
anisotropy parameter f3 z defined as /3 Z = 1 — v^/v^ is 
zero. The models are thus limited to address the degen¬ 
eracy betwe en velocity anisot ropy and dark-halo shape 
as shown by iCappcllarl (120081 ), so relaxing this assump¬ 
tion and considering a non-zero f3 z are both useful to 
enhance the confidence of the conclusions from axisym¬ 
metric mass models. Thus, in order to obtain a more 
realistic and reliable mass distribution of the dark halos 
in dSphs, we employ the axisymmetric mass models that 
are more generalized than those in HC12 and apply the 
models to seven MW and five M31 dSphs. 

This paper is organized as follows. In §2, we briefly de¬ 
scribe the photometric and spectroscopic data used for 
our work and the method of data analysis. In §3, we 
explain our axisymmetric models for density profiles of 
stellar and dark halo components based on a Jeans anal¬ 
ysis and the method of a maximum likelihood analysis to 
be applied to the data. In §4 we present the results of our 
maximum likelihood analysis. Finally, in §5 we discuss 
our results and the implications for dynamical evolution 
of subhalos and association with star formation activity 
in dSphs, and present our conclusions. 
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2. DATA 

In this section we present the data that is used for 
fitting our dynamical models. Table 1 lists the ob¬ 
served properties of seven MW and five M31 dSph galax¬ 
ies: the number of member stars, the stellar mass of 
a galaxy assuming a stellar mass-to-light ratio ( M/L ) 
of unity, V-band absolute magnitude, distance from the 
Sun, projected half-light radius, projected axial ratio 
of a stellar system, lookback time at achieving 70% 
of current stellar mass in its SFH (as detailed in Sec¬ 
tion 5.2), and their references. For the kinematic data 
of their member stars, we use published data as fol¬ 
lows. For Carina, Forn a x, Scul ptor, and Sextans dSphs, 
we adopt IWalker et a]j_ (2009a.bj, a nd for Draco, Leo I, 
and L eo II we adopt iKlevna et all (|2002f ). IMateo et al'l 
<|2008T) , and iKoch et all ( 20(17bth respectively. For the 
M31 dSphs, we select five dSphs (And I, And II, And III, 
And V, and And VII) such that more than 50 stars which 
can be used as kinematic data, because for a dSph with 
less than 50 sample stars it is difficult to evaluate its 
two-dimensional velocity dispersion maps correctly as de- 
sc ribed below. Fo r And II we use the kinematical data 
in Ho et al.l (120121) , and for the other dSphs in M31, those 
in Tollcrud et al.l (1201211 are adopted. 

The kinematic data that we use here are line-of- 
sight velocities taken from the above cited papers. The 
methodology to remove foreground contamination (i.e. 
the Galactic dwarf stars) and reliably identify member 
stars differs in each paper. For Carina, Fornax, Sculp¬ 
tor, and Sextans dSphs, their member stars are estimated 
by ‘exp ectation-maximization’ method in IWalker et alJ 
f 2009m . For Draco, the separation of the member stars 
from the Galactic contaminant stars is clearly made, 
so there is little likelihood of non-Draco stars being 
included in samples. Finally, for Leo I, Leo II, and 
selected M31 dSphs, the method for evaluating mem¬ 
bership is basically a simple kinematical selection, i.e., 
the member stars of those dSphs are determined f rom 
their v elocity distribution (Leo I: IMateo et al.l 120081 
Leo II: IKoch et all 1200 7al And II: IHo et all 120121 the 
other M31 dSphs: iTollerud et aLll2012t) . which is well dis¬ 
tinguished from that of the Galactic foreground stars. 

In order to estimate the linc-of-sight velocity disper¬ 
sion profiles, we adopt here the standard approach of 
using binning profiles. To begin, since we suppose an ax- 
isymmetric system in this work, we analyze the velocity 
data by folding up the stellar distribution into the first 
quadrant in each dSph. In HC12, they estimated line-of- 
sight velocity dispersion profiles along three axes: major, 
minor, and intermediate axes (the last of which is defined 
at 45° from major axis; see Section 3.2 in HC12 for fur¬ 
ther details). However, their method does not utilize all 
of the available member stars in each dSph asa they focus 
on the data only along these three axes, thus their pro¬ 
files may not fully reflect all of the data that is actually 
available for each galaxy. To avoid this, we utilize all of 
the velocity data and derive two-dimensional distribution 
of line-of-sight velocity dispersions as follows. First, we 
represent stellar distribution in the sky (of the first quad¬ 
rant of each dSph) in two-dimensional polar coordinates 
(r, 9), where 9 = 0 is set along the major axis, and then 
divide this into three areas in increments of 30° in the 
direction from 9 = 0° to 90°. For convenience, we call 


these three areas the major axis area, which is defined in 
the region between 0 = 0° and 30°, the intermediate axis 
area with 9 = 30° — 60°, and the minor axis area with 
9 = 60° — 90°. Second, for each area, we radially separate 
stars into bins so that a nearly equal number of stars is 
contained in each bin: ~ 100 stars/bin for Fornax, ~ 80 
stars/bin for Carina and Sculptor, ~ 50 stars/bin for 
Sextans and And II, ~ 20 stars/bin for Draco, Leo I, 
Leo II, and And VII, and 10 ~ 15 stars/bin for And I, 
And III, and And V. We thus derive the velocity disper¬ 
sion maps by using the velocity data of stars contained 
in each bin. The above literatures have considered the 
effect of binary systems on line-of-sight velocity disper¬ 
sions and concluded that such systems have inflated the 
velocity dispersion in each bin by at most ~ 10 %, which 
is small compared to the typical error in the individual 
velocity dispersion, so that the influence of binary sys¬ 
tems in dSphs is in fact negligible. 

Filled circles in Figure [Q display the velocity disper¬ 
sion profiles in three axis areas for 12 dSph satellites. It 
is found that velocity dispersion profiles show systematic 
changes from the galactic center and some notable differ¬ 
ences between the three areas, so that these information 
can distinguish the difference in the mass models. We 
note that in this figure, each data position along the ab¬ 
scissa corresponds to an average galactocentric distance 
of stars within each bin in each of the three axis areas 
for the purpose of demonstration. But in our analysis, 
we actually evaluate an average (r, 9) position of stars 
within each bin in three axis areas and utilize this aver¬ 
aged value in (r, 9) when we fit our mass models to the 
observational data. 


3. MODELS AND ANALYSIS 


3.1. Jeans equations for axisymmetric systems 

The luminous part of dSphs is not really spherically 
symmetric, nor are the shapes of dark matter halos pre¬ 
dicted by high-resolution ACDM simulations. We as¬ 
sume axisymmetry in both stellar and dark-halo compo¬ 
nents, for which axisymmetric Jeans equations are ap¬ 
plied using the velocity dispersion components of stars, 
( v ri v % v z)i cylindrical coordinates. In contrast to 
HC12, we adopt a non-zero velocity anisotropy parame¬ 
te r, 0 Z = 1 — vl/vf, a nd assume Q z = constant as adopted 
in iCappellaril (120081! for the sake of simplicity. This as¬ 
sumpti on is actually sup porte d by recent IV-body simu¬ 
lations; IVera-Ciro et al.l (120141! showed that ACDM sub¬ 
halos have an almost constant /3 2 along the minor axis, 
and only a weak trend as a function of distance along the 
major axis. In this case, the relation between the dark 
matter halo potential, <f>, and moments of the stellar dis¬ 
tribution function is expressed via the Jeans equations: 


v 


2 


1 

v{R,z) 



5 $ , 
z/—— dz, 


(i) 


V 


2 

0 


1 

1-/3, 



R d{vv1) 
7 dR 


R 


dR ’ 


( 2 ) 


where v describes the three-dimensional stellar density. 
For simplicity, we assume that the density distribution 
of the stellar system has the same orientation and sym¬ 
metry as those of the dark matter halo. These veloc- 
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ity dispersions are provided by the second moments that 
separate into the contribution of ordered and random 
motions, as defined by v 2 = <r 2 + v 2 . Since, these equa¬ 
tions indicate intrinsic quantities while we only have pro¬ 
jected moments for the dSphs, we derive the projected 
velocity dispersion from vl and v^, taking into account 
the inclination of the object wi th respect to the o b server , 
following the method given in iTempel fc Teniesl (12006H . 
In the appendix, we describe the detailed model proper¬ 
ties of line-of-sight velocity dispersion profiles, especially 
regarding the effects of velocity anisotropy on the axial 
ratio of a dark halo, thereby we focus on only the main 
properties of our model results in what follows. 


techni ques with the standard Metropolis-Ha sting algo¬ 
rithm ([Metropolis et al.lll953t lHastingslll970H . First we 
define the likelihood function as 

L(M) = IJ^(Oi|M). ( 6 ) 

i 

Here P,(Oi|M) is the conditional probability of finding 
observables CL given a set of model parameters M. Cb 
represents the observables for the ith source, a descrip¬ 
tion of which is provided later, while M represents the 
halo model parameters described in the previous section: 
M = (Q, 6haio, Po, Pz, Oi, i). To construct the conditional 
probability Pj(Oj|M), we suppose the conditional prob¬ 
ability distribution as follows, 


3.2. Stellar and Dark Halo Density Models 

Stellar surface densi ties of dSphs a re empirically fit 
by a Plummer profile (lPlummeilll911jh Thus, here we 
assume this profile and generalize it in the following 
axisymmetric form, v(x,y ) = T(7t& 2 ) -1 (1 + m'f/h 2 )~ 2 , 
where m'f = x 2 + y 2 /q' 2 : q' is the projected axial ratio, 
and (x, y) are the coordinates aligned with the major and 
minor axes, respectively. We estimate &* as the projected 
half-light radius. One of the benefits in the Plummer pro¬ 
file is that one can recover the analytic three-dimensional 
density distribution as follows: 


v(R, z) 


3 L 
47r& 3 




( 3 ) 


where m 2 = R 2 + z 2 /q 2 , so v is constant on ellipses 
with an axial ratio q , and L is the total luminosity. The 
intrinsic axial ratio q is related to the projected ratio q' 
and the inclination angle i such as q' 2 = cos 2 i + q 2 sin 2 i, 
where i = 90° when a galaxy is edge-on and i = 0° for 
face-o n. The int rinsic axial ratio can be derived from 
q=Vq iT — cos 2 i/ sin i, and thus the inclination angle is 
limited by q' 2 — cos 2 i > 0. 

For the dark halo in each dSph, we consider the fol¬ 
lowing density profile 


p(R,z) = Po 



a. 




-(<*+3)/2 


( 4 ) 


f (^!os,i &t,i ) 

2 Var(tTi os>i ) 

( 7 ) 

where eri OSi j is a line-of-sight velocity dispersion corre¬ 
sponding to observables Cb, and Var(<Ji OS; j) is defined 
as the square of the variance of a line-of-sight velocity 
dispersion. N is the number of stars with velocity mea¬ 
surements in the dSph with line-of-sight velocity mea¬ 
surements and a 2 b is the theoretical velocity dispersion 
derived from Jeans equations. For fitting, the velocity 
dispersions of our models are evaluated at the average 
two-dimensional position of each bin associated with the 
observed stars. 

In a practical MCMC method, we calculate the like¬ 
lihood L(M) based on the equation ([5]) for the current 
set of model parameters M. Then, the next set of M' 
is calculated by adding small random fluctuations to the 
previous M, and the likelihood L'(M') is calculated. If 
L'(M')/L(M) > 1, then the next parameter M' is ac¬ 
cepted. If not, we draw a random variable U, which has 
a uniform probability between 0 to 1, and we accept M' 
in case of L'(M')/L(M) > U. If L'(M')/L(M) < U, 
the next parameter M' is rejected and the parameter 
set remains at the previous one M. These procedures 
are iterated for a large number of trials (at least ~ 10 5 ) 
because early trials may retain the effects of initial con¬ 
ditions, which is called the initial “burn-in” phase. 


Pi(aio S ,i\M) = 


V / 27rVar(cri OSii 


: exp 


m 2 = R 2 +z 2 /Q 2 : (5) 

where po is a scale density, fohaio is a scale length in 
the spatial distribution, and Q is an axial ratio of a 
dark matter halo. For simplicity and to focus on only 
an inner profile of the dark matter halo, we suppose 
that the density distribution at outer parts is fixed to 
be p cx r~ 3 . The vantage point of this assumed halo 
model is that it is straightforward to calculate the grav¬ 
itational force, en cl osed mass, and rotation velocity (see 
Ivan der Marel et al.lll994l : iBinncv Sz Tremaindl2008l f . 

In this work, in order to be determined by fitting to 
the observed line-of-sight velocity dispersion, we adopt 
six parameters (Q , frhaio, Po, Pz, ol, *) for each dSph. 

3.3. Markov-Chain Monte Carlo Analysis 

Our aim is to investigate the properties of a dark 
matter halo by exploring the best-fitting parameters 
to the kinematical data of the observed dSph. To do 
this, we utilize Markov-chain Monte Carlo (MCMC) 


4. RESULTS 

In this section, we describe the best-fit dark halo model 
for each of seven MW dSphs (Carina, Fornax, Sculp¬ 
tor, Sextans, Draco, Leo I, and Leo II) and five M31 
dSphs (And I, And II, And III, And V, and And VII) 
and present confidence maps between each parameter. 
Then, we investigate the impact of photometric observa¬ 
tional errors and limited sample volume on the best-fit 
parameters. 

4.1. Best-fit models of dark halos 

We perform the MCMC fitting method for the ob¬ 
served map of line-of-sight velocity dispersions in each 
dSph mentioned above. The best-fit results for each dSph 
are summarized in Table 2. Solid lines in Figure [T] show 
the best-fit line-of-sight velocity dispersion profiles along 
major, minor, and intermediate axes. To begin, we inves¬ 
tigate if the stellar distribution is either oblate or prolate 
as judged from the fitting results. Column 8 in Table 2 
shows the ratio of the most likely value for the prolate 
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TABLE 2 

Results of MCMC analysis for 12 dSph galaxies. 


Galaxy 


Q 

fehalo(pc) 

■£> 

0 

§ 

0 

pc 3 ) 

A 

a 

i (deg) 

-^prolate /-^oblate 

MW dSphs 



709.7±4°® 





—0 09” 1 ” 0 09 

u,u -0.05 

87.ll!;! 


Carina 

0.33 

± 0.02 

0.107 

± 

0.006 

-0.05 ±0.06 

0.0056 

Fornax 

0.38 

± 0.03 

991.1^21 4 

0.086 

± 

0.003 

—0 17+ 016 

u,± ' -0.07 

0.00_o.o4 

90.0_io.6 

0.28 

Sculptor 

0.45 

± 0.03 

637.71^6 

0.168 

± 

0.008 

—0 03" 1 " 0 06 
u.uo_o 04 

O.OO-o.og 

87.8lg! 

0.000015 

Sextans 

0.53 

± 0.06 

1126.7+^'® 

0.028 

± 

0.008 

0 23"*" 012 
u-zo -0.18 

0.00_o.io 

89.8l? 2 2 5 

0.54 

Draco 

0.40 

± 0.05 
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Fig. 1.— Best-fit profiles of line-of-sight velocity dispersions along major, minor, and intermediate axes. The upper-left three panels 
are the cases for Carina, Fornax, and Sculptor, while upper-right three are for Sextans, Draco, and Leo I. The lower-left three panels 
are for Leo II, And I, and And II, and the lower-right three panels are for And III, And V, and And VII. Red, green, and blue marks 
denote observed line-of-sight velocity dispersions along major, minor, and intermediate axes, respectively. Red, green, and blue lines are 
the corresponding best-fit model results along the respective axes. 


case to that for the oblate one. It is found that for most 
of the dSphs, the oblate case yields a much better fit 
than the prolate one, whilst the ratios for some dSphs 
indicate the prolate one is be t ter, es pecially for And II 
dSph. Intriguingly, iHo et al.l (120121 1 investigated kine- 
matical properties of And II and found that this dSph is 
a prolate rotate system, implying that the stellar system 
itself is prolate . Thus, our results are consistent with 
I Ho et all (120121) . In what follows, we use the best-fit 
parameters in the better fitted stellar distribution. 


Before describing the main results, we should inspect 
the influence of the use of different bin sizes on our fit¬ 
ting analysis. To do this, we derive the velocity disper¬ 
sion of three different bin sizes by the same manner as 
described in Section 2. For example, we generate the 
dispersion maps for Sculptor dSph using three different 
number densities of stars in each bin: 45, 80, and 110 
stars per bin, respectively. Then we run MCMC analysis 
for these three types of line-of-sight velocity dispersion 
data. Taking the case of Sculptor dSph as an example, 
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we obtain the difference in the best-fit parameters among 
three different bin sizes is only less than 10 %. In partic¬ 
ular, the best-fit shapes of the dark halos appear to take 
almost the same values: Q = 0.48, 0.45, and 0.47 for 
these three bin sizes, respectively. This is explained by 
the fact that the velocity dispersion profiles are basically 
the same among these different types of data sets. Thus, 
we consider that the difference in the adopted bin size 
has little effect on the results of our fitting analysis. 

Comparing with the previous results of HC12, which 
studied the cases for Carina, Fornax, Sculptor, Sextans, 
Draco, and Leo I, the values of the best-fitting halo pa¬ 
rameters in this work remain roughly the same for these 
galaxies. While we expect a finite degeneracy between 
Q and (3 Z in such a manner that the effect of increasing 
(3 Z in the fitting of velocity disper sions is similar to that 
of decreasing Q f Cappcllarili2008l f. our results as shown 
in Table 2 indicate that most plausible cases for these 
dSphs yield /3 Z ~ 0, so that a flattened dark halo charac¬ 
terized by a small Q is preferred. Thus, our conclusion 
of non-spherical dark halos as the best-fit cases remains 
unaltered. Also, it is clear from Column 2 in Table 2 
that the shapes of dark halos in the MW and M31 dSphs 
can be either oblate with Q < 1 or prolate with Q > 1. 
This result for a prolate dark halo can be understood as 
follows, whereas the case for an oblate dark halo has al¬ 
ready been explained in Section 4.1 of HC12. As shown 
in FigureQ] for dSphs with a prolate dark halo like And I 
and And V, the line-of-sight velocity dispersion along the 
minor axis of these galaxies is relatively small compared 
with those along the other axes. In order to reproduce 
this feature, the best-fit model should yield Q > 1 so that 
v\ decreases in the inner parts. This is because a larger 
Q yields a weaker gravitational force in the z-direction, 
thereby giving decreasing v? in inner parts, as suggested 
from the form of equation [T] Also, the cri os profile along 
the minor axis largely reflects u? profile, so the observed 
small <Ti os along the minor axis suggests Q > 1. We how¬ 
ever note that because of the limitation of the available 
data volume, Qs for these galaxies may be greatly over¬ 
estimated, as we discuss later. 

We investigate other degeneracies in the model 
fitting for determining these six parameters 
(Q, &haio! Po, Pz, ot, i), especially the relation between 
/ 3 Z and the other parameters. In Figure [2] and [3j we 
present 68 % (ler), 95 %, and 99 % confidence levels 
of contours in the two-dimensional plane of Q — /3 Z , 
Po - Pz, a- - p z , fehaio - Pz, and * - f3 z for all sample 
dSphs: Carina, Fornax, Sculptor, Sextans, Draco, and 
Leo I are shown in Figure [21 and Leo II, And I, And II, 
And III, And V, and And VII are shown in Figure [3] 
“FORBIDDEN REGION” indicates the region where 
line-of-sight velocity dispersions are unphysical, namely 
the square of <ti os has negative values in extreme cases 
such as a very flattened dark halo or very large p z . From 
these figures, p z appears to have little degeneracies 
with the other parameters. In particular, the concerned 
degeneracy between Q and p z for the bright dSphs 
having the large number of member stars turns out to 
be rather weak. Whilst this degeneracy certainly exists 
for line-of-sight velocity dispersion profiles along the 
major axis, as described in the appendix, it is broken 
along the minor axis; velocity dispersion profiles along 


this axis are actually more sensitive to Q than /3 Z . This 
break of the degeneracy along the minor axis allows us 
to determine Q, provided the sufficient number of stars 
is available along this axis, as guaranteed for bright 
dSphs. In this respect, it is worth noting that for most 
of Andromeda’s satellites, as well as Leo II, which is 
very faint, and Sextans, for which global luminosity 
distribution is yet uncertain because of its very large 
apparent size, it is impossible to obtain a convergence 
in p z due to a lack of data sample. Similarly, for an 
inclination angle, i, the data in the sample dSphs can 
be reproduced in a wide parameter range of i, since the 
change of i has little influence on line-of-sight velocity 
dispersion profiles compared to the change of the other 
parameters. Thus, with available data alone it is difficult 
to determine the best-fit values for both p z and i. 

From our best-fit results, the inclination angles for all 
dSphs are over 70°, that is, all stellar systems are nearly 
edge-on. One of the possible reasons for this result is 
the limitation of an inclination angle from the relation, 
q 12 — cos 2 * > 0. This inequality suggests that for a given 
projected axial ratio, q ', an inclination angle is restricted 
to be larger than cos -1 q'. It is also worth noting that 
in M31’s dSphs for which only a small number of spec¬ 
troscopic data is available, the confidential intervals of i 
cover rather a large part of parameter space, in contrast 
to the MW dSphs being confined in high inclination an¬ 
gles. 

4.2. Central density profiles of dark halos 

In the following, we focus on central density pro¬ 
files of dark matter halos for the sample dSphs. The 
question of which central density profiles, cusped or 
cored, observed dSphs have has been ambiguous, and 
thus lau nched the debate known as a core-cusp prob¬ 
lem. iGilmore et al.l (2007 ) gave a review of studies on 
the inner slope of dSphs’ dark halo profile based on 
spherical Jeans modeling and concluded that cored pro¬ 
files are preferred for all classical dSphs in the MW. Re¬ 
cently, applying the relation betw een the ha lf-lig ht radii 
and m asses within these radii, IWalker fc Penarrubial 
12011 1 inferred the central density profiles in Fornax 
and Sculptor. They suggested that Navarro-Frenk- 
White (NFW)-like cuspy profiles can be ruled out at 
significance levels > 96% and > 9 9% for Fornax 
and S culp tor, respectively. More over, lAmorisco et al.1 
(l20l! and lAgnello fc Evansl (|2012l f applied the projected 
virial theorem to the multiple stellar components for For¬ 
nax and Sculptor and argued that this modeling disfa¬ 
vors cuspy profiles in these dSphs. I n contrast, us ing the 
axisynm ietric Schwarzsc hild method. I Jardel et al.l ([2013H 
and i.Tardcl fc Gebhardtl (1201311 applied their axisymmet- 
ric stellar models to MW dSphs. They concluded that 
the dark halo profiles in the dSphs is similar to NFW 
profiles. iStrigari et al.l (12014H assumed the distribution 
function of the stellar system embedded in a spheri¬ 
cal NFW dark halo, and showed that this distribution 
function reproduces the line-of-sight velocity dispersions 
and surface stellar densities of two subcom ponents i n 
the Sc ulptor dSph. On the o t her h and, IWalker et al.1 
(l2009dl and iBreddels fc Helmil (1201311 concluded it dif¬ 
ficult to constrain an inner profile of a dark halo with 
current data quantity and quality. 

We point out that all of the previous studies have as- 
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TABLE 3 

Comparison of Best-fit Parameters Estimated from Full Data Sample with those Estimated from 

Data Only within Half-light radius 


Object 

Data 

Q 

bhalo(pc) 

Po (M 0 pc ®) 

P. 

a 

i (deg) 

Fornax 

Data within r^aif 
Full data 

0.83 ± 0.04 

0.38 ±0.03 

708.5!??' 7 

OOLltaliS 

0.101 ±0.004 

0.086 ± 0.003 

0.40 ± 0.03 

—0 17+ 016 
-0.07 

-0.14 ±0.06 

0 . 00 _o.o 4 

89'5l? 0 5 6 

90.0_io.6 

Sculptor 

Data within 7 + a if 

O 

^1 

CO 

1 + 

o © 

o o 

450 . 71 ®).;® 

0.241 ±0.014 

0.48 ± 0.05 

—0 02”*” 0 ' 02 
u.u ^_ 0>06 

89.5t% s 6 


Full data 

0.45 ±0.03 

637 7+ 32,6 

vot . ( _26.0 

0.168 ±0.008 

—0 03 +0 ' 06 
U.UO-o 04 

0 . 00 _o.o 9 

87.81®;? 


TABLE 4 

Best-fit Parameters for Considering Photometric Measurement Errors 


Object 

Measurement Error 

Q 

£>halo(pc) 

P0 (Mg pc ®) 

/3a 

a 

i (deg) 

Fornax 

Half-light radius 








Upper limit 

0.39 ± 0.02 

941.31?® g 

0.087 ±0.005 

-0.13 ± 0.08 

0 . 00 _o.o 4 

90.0_io.5 


Lower limit 

0.40 ± 0.02 

924.71?®;® 

0.094 ±0.006 

-0.10 ±0.09 

-0.04 ±0.04 

88 - 61 );) 


Projected axial ratio 



-o-iel?;)? 

—0 01+ 0 01 
U-U-L-o.03 


Upper limit 

0.40 ± 0.02 

887.91?®;) 

0.095 ± 0.004 

90.0—12.1 


Lower limit 

0.37 ± 0.02 

1013.91??'® 

0.088 ± 0.004 

-0.09!°;“ 

0 . 00 _o.o 5 

90.0-10.8 

Sextans 

Half-light radius 








Upper limit 

0.69 ± 0.08 

1110 . 61 ?®'? 

0.024 ± 0.002 

0.37 ±0.11 

0 . 00 _o.o 9 

90.0_i2.o 


Lower limit 

0.36 ±0.05 

1410.ll?®?;? 

0.026 ± 0.003 

< 0.16 

—0 in + 9 - 98 
u - lu -0.06 

89.7l? 2 ®i 


Projected axial ratio 






Upper limit 

0.78 ± 0.08 

896.71®,?') 

0.030 ± 0.003 

0.38!?;?) 

—0 io +01 ° 

u.iu_ 0 08 

89-8l? 4 2 8 


Lower limit 

0.48 ± 0.06 

1047.71?)'® 

0.033 ± 0.003 

0 29" 1 " 0 ' 11 
u ' zy -0.13 

—0 02”*” 0 02 
u ' uz -0.08 

89.51?+ 

And II 

Distance 

Upper limit 

0.51 ± 0.03 

2886.ll?®®'® 

0.0085 ± 0.002 

0.17 ±0.08 

—0 01 + 0 - 01 
u - ul -0.07 

89-81?+ 


Lower limit 

0.53 ± 0.03 

2818.41?)?;® 

0.0093 ± 0.002 

0-28!?;?) 

0.00-0.042 

89.7l? 2 ® 5 

And VII 

Distance 





Upper limit 

-| 07 +O .45 
-0.30 

4ii-5l??;? 

0.078 ±0.012 

-0.03!?;®® 

—0 76~*”°‘ 27 

u - ‘°-0.31 

70.0l??'? 


Lower limit 

1 61+ 0 - 38 
i - D± — 0.27 

480.51®?;® 

0.068 ± 0.01 

-0.051?;®® 

—0 61 +0 ' 24 
U-D1 — 0.30 

70 i + 19 - 9 
* — 14.6 


sumed a spherical dark halo. Thus, our non-spherical 
dark halo models are expected to set useful constraints 
on central profiles of dark halos. From our fitting re¬ 
sults, we find that not all of the dark halos in the sample 
dSphs have a cored central density profile; while most 
of the dSphs indicate constant density profiles or shal¬ 
lower cusps, the Draco, Leo I, And II, And III, And V, 
and And VII dSphs show a steep inner density slope, 
a, with -0.861°;??, -1.40±g;g|, -1.03™, -1.43!°;?), 

—1.33!?;??, and — 1.341?;??, respectively. Therefore, in 
the light of axisymmetric mass models combined with 
currently available data set, these dSphs suggest the pres¬ 
ence of NFW-like or more strongly cusped dark halos. 


4.3. The impact of a sample selection and observational 
errors on best-fit parameters 

As indicated above, we find that dSphs in the MW 
and M31 have non-spherical dark halos, where large and 
diffuse dSphs are characterized by low dark matter den¬ 
sities. However, it is unclear how the ways of data sam¬ 
pling (such as data volume and distribution of member 
stars) as well as observational errors (such as in the half- 
light radius, projected axial ratio of luminous parts, and 
distance from an observer) affect the estimation of the 
halo parameters. Therefore, in what follows, we take 
into account these effects in the likelihood analysis and 
investigate their impact in the best-fit dark halo param¬ 


eters. 

First, in order to inspect the effects of data volume 
and distribution of member stars, we derive line-of-sight 
velocity dispersions using the data only within half- 
light radii for Fornax and Sculptor dSphs, and then run 
MCMC analysis. Since these galaxies have the largest 
number of data sets and their member stars are widely 
distributed beyond their half-light radii, these are suited 
to assess the effects of the limited data volume and dis¬ 
tribution of member stars. As shown in Table 3 and 
Figure O we find that the best-fit halo parameters using 
the data only within the half-light radii are substantially 
different from those using the full data sample. It is also 
found that using the data only within the half-light radii, 
the degeneracy between Q and /3 Z has emerged signifi¬ 
cantly compared to using the full data sample. There¬ 
fore it is difficult to determine these parameters due to 
the lack of data sample in the outer regions. This is be¬ 
cause, as described in the appendix, the impacts of /3 Z on 
the line-of-sight velocity dispersion along the major axis 
appear prominently at outer parts. Therefore, if a suf¬ 
ficiently large number of sample stars is available in the 
outer region of a dSph, we can obtain tight limits on (f z 
with small uncertainties. Consequently, the best-fit pa¬ 
rameters for the dSphs in which spectroscopic informa¬ 
tion are incomplete in their outer region (e.g., Sextans, 
And I, III, V, and VII) are subject to the effects of data 
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Draco Leo I 



Fig. 2. — Likelihood contours for each dark halo parameters against velocity anisotropy, (3 Z , for six dSphs. Clockwise from top left, 
Carina, Fornax, Sextans, Leo I, Draco, and Sculptor. Contours show 68% (lcr), 95%, and 99% confidence levels. The cross point in each 
panel indicates the best-fit value of each parameter. Dotted line in each panel is a boundary with “FORBIDDEN REGION” in which 
line-of-sight velocity dispersions are unphysical (see the text in Section 4.1 for details). 


deficiency, and thus we suggest that in order to derive 
more confidential dark halo structure in these dSphs, we 
require the observational data over much larger areas. 

Second, we repeat the MCMC analysis to obtain the 
best-fit parameters considering uncertainties in the half- 
light radius, rhaif, and projected axial ratio, q ', of the 
stellar system, because what stellar density profile is 
adopted affects the best-fit parameters of a dark halo 
through Jeans equations (lEvans, An fe Walkerl2009l) . To 
investigate this, we estimate best-fit parameters for For¬ 
nax and Sextans by adopting the upper and lower obser¬ 
vational limits of rhaif and q' within their uncertainties. 
These dSphs have similar rhaif and q' but their errors 
in Sextans are larger than those in Fornax as shown in 
Table 1. Table 4 and Figure [5] show the effects of the er¬ 
rors in rhaif on dark halo parameters. For Fornax, with a 
smaller standard error in rhaif, fitting results are hardly 
affected by this error, while a relatively large error in 
rhaif for Sextans affects the best-fit Q and p z . On the 


other hand, the effects of observational errors in q' are 
weak compared with those of rhaif as shown in Figure [G] 
We note that for the t>haio — Pz and po — Pz relations in 
Fornax, each contour deviates rather largely because of 
the degeneracy between po and 6haio- Likewise, we inves¬ 
tigate the impact of measurement errors in the distance 
from the Sun using the data of And II and VII dSphs 
which show large uncertainties in the distance estima¬ 
tion compared with MW ones (see Table 1). As shown 
in Figure [3 we find that the distance error as reported 
for these dSphs does not affect the results of maximum 
likelihood analysis. 

Accordingly, the constraints on dark halo structures in 
dSphs are affected largely by the lack of kinematic sample 
and distribution of member stars rather than uncertain¬ 
ties in photometric data. Thus, to set robust constraints 
on dark halo structures in dSphs, we require deep photo¬ 
metric data to assemble many sample stars down to faint 
magnitudes and spectroscopic data over large areas out 
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Fig. 3. — The same as Figure[2]but for the other dSphs. Clockwise from top left, Leo II, And I, And III, And VII, And V, and And II. 


to the tidal radii of dSphs. 

5. DISCUSSION AND CONCLUDING REMARKS 
5.1. Comparison with A CDM subhalos 

Recently, lYera-Ciro et al.1 (12014 hereafter VC14) in¬ 
vestigated the distribution for axial ratios of triaxial sub¬ 
halos obtained from the Aquarius simulations. They 
found that axial ratios of subhalos with maximum cir¬ 
cular velocity in the range 8 < V max < 200 km s -1 are 
little different between isolated subhalos and those asso¬ 
ciated with host galaxies. Moreover, these subhalos can 
be approximated as oblate axisymmetric objects; the av¬ 
erages of intermediate-to-major axial ratios, ( b/a ), and 
minor-to-major axial ratios, (c/a), (supposing a > b > c) 
at r ~ 1 kpc are (6/a) ~ 0.75 and (c/a) ~ 0.60, respec¬ 
tively. They also analyzed the velocity structure in the 
cylindrical radial and vertical directions, i.e., ctr and a z , 
respectively, estimated anisotropy /3 Z = 1 — cr^/cr^, and 
computed this as a function of distance along minor and 
major axes. Their analysis shows that along the mi¬ 
nor axis a z is almost the same as ctr, that is, (3 Z ~ 0, 


whereas along the major axis /3 Z > 0 near the halo cen¬ 
ter and (3 Z < 0 at the outer parts. This systematic trend 
is, however, weak with large scatters among the individ¬ 
ual subhalos, thus, the radial dependence of 0 Z may be 
regarded relatively weak. 

Although these simulation results may support the va¬ 
lidity of our axisymmetric models for the dark halos of 
dSphs, the axial ratios of dark halos in both MW and 
M31 dSphs are systematically smaller, and thus more 
flattened, than those of ACDM subhalos. We note that 
for this comparison with VC14’s simulations, in the case 
of prolate halos (Q > 1), where the minor axis of a dark 
halo is aligned with the major axis of a stellar distribu¬ 
tion, we use Q' = 1/Q instead of Q following VC14, so 
Q' = 0.41 for And I and 0.21 for And V. As shown in 
Figure 7 in VC14, there are no subhalos with axial ratios 
less than 0.6, which accords with the interpretation that 
less-massive halos tend to be more relax ed and there¬ 
fore m ore s pherical than massive halos. IKuhlen et al.l 
(I2007f) and ISchneider et al.1 (12012H also have concluded 
from their TV-body simulations that less-massive dark ha- 
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Fig. 4. — Comparison of likelihood contours calculated from dif¬ 
ferent data selections for Fornax (upper panel) and Sculptor (lower 
panel). Contours drawn with solid and dashed lines show 68% (lcr), 
95%, and 99% confidence levels calculated from full data sample 
and data only within their half-light radius, respectively. 




Fig. 5.— Comparison with likelihood contours considering mea¬ 
surement errors in half-light radius for Foranx (upper panel) and 
Sextans (lower panel). Red, green, and blue contours show 68% 
(lcr) and 95% confidence levels for the best-fit case (as shown in 
Figure EJ, the cases considering an upper and lower limit in the 
error of r^aif* respectively. 

los have mildly triaxial shapes. In contrast, our present 


Fig. 6.— Same as Figure [5] but for considering measurement 
errors in the projected axial ratio of the luminous component. 




Fig. 7.— Same as Figure [5] but considering measurement errors 
in the distance from the Sun for And II (upper panel) and And VII 
(lower panel). 

analysis for the dark halos of dSphs suggests that Q is 
much smaller than the prediction of these simulations. 
Thus, there exists a mismatch in the shapes of subhalos 
predicted from ACDM theory. 

One of the possible mechanisms to alleviate this dis- 
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Lookback time (Gyr) 


Fig. 8. — Cumulative star fo rmation history of dwarf satellites in 
the MW and M31 taken from lWeisz et all (120141) . 



%7 (Gyr) 


Fig. 9. — Relation between p(fehalo) and the lookback time achiev¬ 
ing 70% of current stellar masses. 


crepancy may be tidal effects from a host dark halo. 
Recent JV-body simulations show that subhalos be¬ 
come significa ntly elongated at pericenter passage (e.g., 
iKuhlen et al.l 120071) because of undergoing strong tidal 
force from de e p pote ntial of its host dark halo. Moreover, 
iBarber et al.l (120151) predicted that dSphs with lower 
M/Ls have more spherical dark halos due to tidal strip¬ 
ping, whilst more dark matter dominated systems are 
more triaxial. In light of these predictions, it is inter¬ 
esting to remark that the dark halos of Fornax, Sculp¬ 
tor, Leo I, and Leo II with relatively low M/Ls tend 
to be less flattened than Carina, Sextans, and Draco 
with high M/Ls. To quantify this property, we calcu¬ 
late the average of dark halo axial ratios, < Q >, for 
these low and high M/L dSphs, respectively, and obtain 
< Q >iow = 0.65 ± 0.08 and < Q >Mgh= 0.42 ± 0.04. 
T herefore, the trend fo r the shapes of subhalos predicted 
bv lBarber et alJ (|2015l) are roughly consistent with those 
for the shapes of dark halos associated with MW dSphs. 

However, the axial ratio of these simulated subha¬ 
los at the closest pericentric distance is never smaller 


than 0.5, and thus we should consider additional mech¬ 
anisms to make a dark halo more flattened. We sug¬ 
gest that baryons may potentially invoke reproduction 
of very flattened dark halos. First, owin g to adiabatic 
contr action during disk formation (e.g., Dutton ct al.1 
l2007f) . a dark halo with a massive disk galaxy can have 
a deep and steep gravitational potential in its central 
part compared with a pure dark halo without baryons. 
Then, when subhalos pass through this deep potential 
of the host halo, the shape of subhalos may be more 
flattened by strong tidal distortion than the case with¬ 
out considering baryo nic effects in the host halo. Sec¬ 
ond, as envisaged by iBrvan et al.l (120131). while cooled 
baryo ns tend to round a dark halo (jKazantzidis et akl 
I2004D . their removal due to strong feedback in the 
central part s of dSphs, as seen i n recent simulation 
studies (e.g..lGovernato et al.ll201(i iTevssi er et al.ll2013l 


cyssicr c 

Ul2014D, 


iDi Cintio et al.l l2014aU bl: iMadau et all I2U14I ). may pre¬ 

vent dark halos of such dSphs from being rounder, i.e. 
their shapes remain flattened. It is, however, still un¬ 
clear how the combination of these baryonic effects in¬ 
deed modifies the shapes of dark subhalos. To resolve 
these issues, while more observational data are to be as¬ 
sembled for more robust determination of dark halo pa¬ 
rameters, further simulation studies fully taking into ac¬ 
count baryonic effects are worth exploring to get impor¬ 
tant insights into the issue of dark halo shapes addressed 
here. 


5.2. The relation between dark-halo structure and SFH 

As mentioned in the previous subsection, recent nu¬ 
merical simulations imply that internal dark matter 
structure in the low-mass halos can be altered by stellar 
feedbacks associated with star formation activity. There¬ 
fore, it is natural to expect that dark halo properties 
depend on the SFH of the stellar components, and we 
thus investigate whether this dependence indeed exists 
by comparing it with the observed SFH of dSphs. For 
this p urpose, we adopt the recent work by IWeisz et al.l 
(l2014h . which derive the SFH of dwarf galaxies in the Lo¬ 
cal Group based on the analysis of their color-magnitude 
diagram taken from deep imaging of the Hubble Space 
Telescope. Figure[5]displays the cumulative SFH of dwarf 
satellites in the MW (except for Sextans) and M31 taken 
from their work, indicating various formation histories 
such as rapid and consecutive star formation. To facili¬ 
tate the comparison with the present work, we estimate 
the lookback time at achieving 70% of the current stellar 
mass of these dSphs, To .7 (as indicated as a black horizon¬ 
tal line in Figure [ 8 ] and this value in each dSph is shown 
in Table 1); To .7 may characterize the duration and effi¬ 
ciency of star formation in dSphs. This is because if we 
choose 50% instead of 70% of the current stellar mass 
as an indicator, we cannot completely trace the all ac¬ 
tivity phases of star formation that would be capable of 
changing dark-halo structures. For instance, as is shown 
in Figure [ 8 l Fornax, Leo I, and And II have the periods 
of active star formation after the stellar masses exceed a 
half of these current ones. We compare between r 0 .7 and 
dark halo parameters obtained from our analysis. Fig¬ 
ure E shows the relation between p{b h a i 0 ) and To. 7 . We 
consider MW and M31 dSphs separately because they 
are located in the different host halos and there is actu¬ 
ally a systematic difference in quenching times of their 
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sta r formations due to environmental effects as suggested 
by IWeisz ct al.1 (120151) . It is found that dSphs showing 
relatively consecutive SFH, namely lower values of To .7 
(e.g., Leo I, Fornax, and And II) have a low dark matter 
density within a central region, whereas those with rapid 
SFH (e.g., Sculptor, Draco, and And VII) have a dense 
and concentrated dark halo. We interpret that the cen¬ 
tral density of a dark halo may be gradually decreased 
associated with gaseous outflow driven by effects of stel¬ 
lar feedbacks. 

We also inspect other possible relations between r 0 .7 
and dark halo parameters such as an axial ratio, Q , and 
a central density slope, a, but there are no remarkable 
relationships within uncertainties of these parameters. 
This is because, as mentioned above, the shapes of dark 
halos would be changed by not only stellar feedbacks, 
but also tidal effects from a host galaxy and dark halo. 
On t he other hand, for centra l dens ity slopes of a dark 
halo, iNipoti fe Binnevl (12015T) and iDel Popolo fe Pace] 
(!2015fl proposed the possible mechanism of the core for¬ 
mation that the main process of transforming cusps to 
cores would be dynamical friction of gas clumps. More¬ 
over, it may give only a small impact on p(b halo), because 
this transformation occurs in the innermost regions of a 
dark halo. According to their calculation, this process 
may occur prior to star formation and subsequent feed¬ 
back processes, implying that the relation between a and 
SFH may be made hazy, in agreement with the observed 
properties. Finally, we note that we find no explicit rela¬ 
tion between dark halo parameters and averaged meta l- 
licity, <[Fe/H]>, which is taken from lKirbv et alJ (1201311 . 
This may be due to the fact that <[Fe/H]> in the stars 
reflect several processes such as a gas outflow as well as 
SFH, so the direct relation with background dark-halo 
properties may be weakened. 

Thus, although both the available data and mass mod¬ 
eling for dSphs are yet limited, we suggest the possible 
link between the density of their dark halos and SFH of 
the dSphs. In order to explore this relation in more de¬ 
tail, we require not only further refined modeling such 
as considering triaxial structures of dark halos, but also 
assembling more photometric and kinematic data of the 
dSphs in their larger parts. 


5.3. Constant dark matter surface density within a 
radius of maximum circular velocity 

Some previous studies have noted an intriguing gen¬ 
eral property of da rk h a los i n a various kinds of g alaxies . 
For the studies of iKormendv fe Freemanl (12004L 120141) . 
using 55 spiral galaxy rotation curves and the line-of- 
sight velocity dispersions of a few dSphs, they evalu¬ 
ated the central density po and core radius r c of dark 
halo assuming the non-singular isothermal sphere mod¬ 
els, and they first found that p$r c is almost indepen¬ 
dent o f galaxy luminosity, pgr r oc l ™^ 0 ' 067 Sim¬ 
ilarly. iGentile et ahl (|2009l) . iDonato et al.l (12009D . and 
iSalucci et al.1 ( 20121) also found and confirmed, respec¬ 
tively, the constant dark matter surface density, despite 
of assuming any cored dark matter density profiles such 
as Burkert and pseudo-isothermal models. These works 
have provided us considerable insight for the universal 
properties of dark-halo structure, and thus we investi¬ 
gate whether this constancy exists in our estimated dark 
halos in dSphs. Figure [TO] displays the relation between 
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Fig. 10.— Product of p(fthalo)&* (upper panel) and p(bhalo)&halo 
(lower panel) as a function of B-band absolute magnitude. The 
solid lines and in each panel denote our best-fit line to the data 
using least squares method. 


dark matter surface densities and luminosities of each 
dSph. The upper panel shows the product of dark mat¬ 
ter densities at its scale length, p(&haio), and half-light ra¬ 
dius of a stellar component, p(&haio)(>*, whilst the lower 
panel is p(&haio)^haio- In order to investigate the con¬ 
stancy of above two relations, we employ a least squares 
fitting method to determine the slope of these surface 
densities as a function of luminosity, and we find that 
the slopes of those are p( 6 haio)fr* oc £^ 045±0 ' 010 , and 
p(frhaio)khaio oc 2 ,-°- 08±0011 ; respectively. It is clear that 
p(bha\o)b * and p^haio^haio are nearly constant values 
with respect to the luminosity of these galaxies. In com¬ 
parison with above previous results, our results have a 
similar constancy, even though assumed dark halos have 
cusped profiles. 

However, because of the difference of assumed dark 
halo density profiles among this comparison, the defi¬ 
nitions of dark matter surface densities of each study 
should also be different. Thus, in order to inspect 
more universal dark halo properties without the am¬ 
biguously of the assumed dark halo density profile, we 
introduce the new definition for dark halo surface den¬ 
sity, namely, the mean dark halo surface density within 
a radius of maximum circular velocity, Sy max . Since 
this surface density was already defined in our recent 
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Fig. 11. — Mean surface dark halo density, Ey max as a function of B-band absolute magnitude. The rhombic points denote the dSphs, 
which ha ve cored (magenta) a nd cusped (cyan) dark halos, based on this work. On the ot her hand, the other s ample data ar e the original 

S ic of Donato ct al. (2009). Blue and red triangles are the nearby spirals in THINGS (dc Blok ct al. 2008) and GHASP dSpano et al.l 
data, respectively. Green circles and black squares are dwarf irregular galaxies ([Gentile et a l. 2005, 2007), and spiral and elliptical 
galaxies investigating by weak lensing (Hoekstra et al. 2005), respectively. The dark halo density distribution of these sample galaxies are 
assumed by any cored profiles such as Burkert and pseudo-isothermal profile. 


paper (jHavashi fo Chibal I2015L hereafter HC15), we de¬ 
scribe only the important properties of this. This sur¬ 
face density is directly correlated with the central dark 
matter surface density, Ey max oc p»r*, where p* and 
r* are central density and scale length of arbitrary 
dark halo profiles, respectively. In HC15, we calcu¬ 
lated Ey max using the data of dark halo profiles based on 
the HI gas rotation curve of late- a nd early -t ype spirals 
with pseudo-isoth ermal dark halos (Ide Blok et al.1120081 : 
iSpano et all I200 8T) . dwarf irreg ulars with Burkert dark 
halos (iGentile et akll2005L I2007T) . the galaxy-galaxy weak 
lensing from spiral and elli ptical g al axies with Burk¬ 
ert profiles (analyzed bv iDonato et al.l [20091 data from 
iHoekstra et al.l I2005T ) and the dSphs in this work. Then 
we found that Ey- max does not depend on the maximum 
circular velocity corresponding with dark-halo mass of 
these galaxies. Figure [TT] shows the Em, half — Mb rela¬ 
tion for same sample in HC15 with luminosity over al¬ 
most 14 orders of magnitude. As shown in this figure, 
we confirm that this dark-halo surface density has also 
constancy against galaxy luminosity without concerning 
the assumption of dark halo profiles. 

5.4. Concluding Remarks 

In this paper, we revisit the dynamical analysis of non- 
spherical dark halos for the dSphs in the MW and M31 
based on revised axisymmetric mass models. In contrast 
to HC12, which assumed h? = v ^ for simplicity; we relax 
this constraint and set a non-zero velocity anisotropy, 


f3 z = 1 — Vz/Vft, as this anisotropy can be degenerate 
with the effect of Q. Based on the application of our 
models to seven MW dSphs and five M31 dSphs, we find 
that the best-fitting cases for most of the dSphs yield 
not spherical, but flattened halos even considering f3 z . It 
is worth noting that the best-fit parameters, especially 
for the shapes of dark halos and velocity anisotropy, are 
susceptible to both the availability of kinematic data in 
the outer regions of the system as well as the effect of the 
small number of sample stars. Thus, to estimate more 
plausible dark-halo structure in the dSphs, we require the 
photometric and kinematic data over much larger areas. 
We have found from the revised analysis and currently 
available data that the derived shapes of dark halos in 
dSphs remain systematically more flattened than those 
of dark matter subhalos calculated from ACDM-based 
IV-body simulations, as suggested in HC12. This mis¬ 
match needs to be solved by a theory including baryon 
components and the associated feedback to dark halos 
as well as by further observational limits in larger ar¬ 
eas of dSphs. It was also found that more diffuse dark 
halos may have undergone consecutive SFH, as charac¬ 
terized by low ro. 7 , implying that the formation process 
of dSphs is affected and thus imprinted in the structure 
of their dark halos. In the near future, planned surveys of 
the MW and M31 dSphs using Hyper Suprime-Cam and 
Prime Focus Spectrograph attached on the Subaru Tele¬ 
scope (jTakada et al.ll2014() will enable us to hunt a num¬ 
ber of faint stars in the outer parts of dSphs and measure 
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their kinematic data and metallicities, thereby allowing 
us to obtain severer limits on the dark halo distribution 
and characterize the dynamical evolution of dSphs. 
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APPENDIX 

EFFECTS OF Q AND [i z ON LINE-OF-SIGHT VELOCITY DISPERSION 

In this appendix, we demonstrate the impact of the non-spherical shape of a dark hal o, Q, and th e veloc ity anisotropy 
f3 z of member stars on their line-of-sight velocity dispersion profiles, motivated by the lCannellaril (120081 1 paper, which 
showed t hat t he variations of Q and 8 Z are a similar effect on cri os profiles. 

Figure [A. II indicates normalized line-of-sight velocity dispersions along the major and minor axes, respectively. We 
set the fixed axial ratios of stellar distribution in this t est calculation, q = 0.7, which is a typical value in MW dSphs. 
First, as is shown in the upper panels of Figure IA.1I denoting the cri os profiles along the major axis, the effects of 
increasing ft z (top left) and of decreasing Q (top right) are resemblant in the property that these weaken the wavy 
feature of the cri os profile caused by a non-spherical Q along the major axis, as already discussed in HC12. However, 
there is a crucial difference between these effects. The difference in crios profiles by the change of @ z does occur at 
not the inner but outer parts. This is because v ^ including a constant velocity anisotropy in Equation (2) contributes 

to i7i os more than dose in outer parts, thus the impacts of fi z on cri os emerges at only outer parts. On the other 
hand, for the case of the change of Q, this is opposed to that of /T. This is because a larger (smaller) Q yields weaker 
(stronger) gravitational force in the z-direction, thereby decreasing (increasing) in inner parts. 

Second, as shown in the lower panels of Figure lATl there is no difference by changing /3 2 , because cri os along the 
minor axis is not contributed by v^, that is, we need not to consider the effects of /? 2 . On the other hand, the effects of 
changing Q are monotonous, which can be straightforwardly understood. As mentioned above, we can consider only 
the impact on v%, thus a decreasing (increasing) Q reduces (increases) vl in inner parts. 

In summary, a finite degeneracy between Q and /3 2 can be broken by investigating a systematic difference in the 
effects of these parameters on tri os profiles between the minor and major axes. This suggests that to definitely determine 
Q , an axial ratio of a dark halo, we need a sufficiently large number of kinematic data for member stars, whereby 
accurate crios profiles along both minor and major axes are available. 
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Fig. A.l.— Upper panels show normalized line-of-sight velocity dispersions, rri/(Gpob*) 1 /2 , along the major axis for a cored case 
(a = 0), whereas the bottom panels show these velocity dispersions along the minor axis for the same case. For all of these cases we 
suppose that the edge-on galaxy (i = 90°) and the ratio of hhalo/£>* i s unity for the sake of demonstration. 


























